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Abstract 



Poles of solutions to the Painleve-I equations are intimately related 
to the theory of the cubic anharmonic oscillator. In particular, poles of 
integrale tritronquee are in bijection with cubic oscillators that admit 
the simultaneous solutions of two quantization conditions. We analyze 
this pair of quantization conditions by means of a suitable version of 
the complex WKB method. 

1 Introduction 

The aim of the present paper is to study the distribution of poles of the 
solutions y = y{z) to the Painleve first equation (P-I) 

y" = Qy'' - z, z&C , 

with particular attention to the poles of the integrale tritronquee. 

As it is well-known, any local solution of P-I extends to a global meromor- 
phic function y{z),z G C, with an essential singularity at infinity |GLSOO| . 
Global solutions of P-I are called Painleve-I transcendents, since they cannot 
be expressed via elementary functions or classical special functions |Inc56| . 
The integrale tritronquee is a special P-I transcendent, which was discov- 
ered by Boutroux in his classical paper |Boul3| (see |JK88| and |Kit94| for 
a modern review). Boutroux characterized the integrale tritronquee as the 
unique solution of P-I with the following asymptotic behaviour at infinity 

/z Att 
-, if I argz| < — . 

Nowadays Painleve first equation is studied in many areas of mathemat- 
ics and physics. Indeed, it is remarkable that special solutions of P-I describe 
semiclassical asymptotics of a wealth of different problems (see |Kap04| and 
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references therein). In particular, in |DGK09| it is discovered that the inte- 
grate tritronquee provides the universal correction to the dispersionless limit 
of solutions to the focusing nonlinear Schrodinger equation. 

Theoretical and numerical evidences led the authors of |DGK09] to the 
following inspiring 

Conjecture. // a G C is a pole of the integrale tritronquee then |arga| > 

5 • 

Following the isomonodromic approach to P-I |Kap04| , any solution y{z) 
gives rise to an isomonodromic deformation of the following linear equation 
with an irregular singularity 



The deformation of the equation is manifestly singular at every pole a G C 
of y, however in Theorem |2] we show that at the singularity this equation 
can be replaced with a simpler one, which has the same monodromy data 
(cf. |IN86| for Painleve II). This is the following Schrodinger equation with 
cubic potential 



Here a is the location of the pole of y and 6 is a complex number entering 
into the Laurent expansion of y around a (see formula ([7| below). 

The isomonodromy property implies that there exists a natural injective 
map Ai from the space of solutions of P-I to the space of monodromy data 
of the above equations (see Lemma [2]) , while the Schrodinger equation de- 
fines naturally a map T from the space of cubic potentials to the space of 
monodromy data. 

Our first main result is Theorem |3l which states that a S C is a pole of 
y{z) if and only if there exists b £ C such that M{y) = 'T{V{X; a, b)). 

In particular, because of the special monodromy data related to the inte- 
grale tritronquee (see Theorem [H due to Kapaev), we will show that the poles 
of the integrale tritronquee are in bijection with the simultaneous solutions 
of two different quantization conditions. 

The above approach naturally embeds the study of poles of Painleve- 
I transcendents into the Nevanlinna's theory of branched coverings of the 
sphere and the complex WKB method. 

The beautiful theory of R. Nevanlinna (see |Nev70| and |Elf34| ) relates 
bijectively the Schrodinger equations with a polynomial potential to the 
branched coverings of the sphere with logarithmic branch points, considered 
up to conformal equivalence. Using this theory we are able to prove the 
surjectivity of the map A4 (see Theorem |5)). 




¥(A) 
dA2 



V{X; a, b)i){X) , V{X; a, b) = AX^ - 2aX - 28b . 
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Moreover Nevanlinna's theory provides the poles of any solution of P-I 
with an unexpected and remarkable rich structure. In particular, poles of the 
tritronquee solution can be labelled by the monodromy of coverings of the 
Riemann sphere with 3 logarithmic branch points. In a subsequent paper, we 
are going to use this topological description to complete the WKB analysis 
of the present paper. 

The WKB analysis of P-I developed in |KT05| has never been applied to 
the direct study of the distributions of poles. To achieve such a goal we fol- 
low the Fedoryuk's approach (see |Fed93| ) to the complex WKB theory, and 
in the Classification Theorem we give a complete topological classification of 
the Stokes complexes of all cubic potentials. As a consequence of the Classi- 
fication Theorem, we obtain our second main result: all polynomials whose 
monodromy data, in the WKB approximation, are the monodromy data of 
the integrale tritronquee have the same topological type of Stokes complex 
and satisfy a pair of Bohr-Sommerfeld quantization conditions, namely sys- 
tem (|25p . In particular, in this way we reproduce the conditions obtained 
by Boutroux, through a completely different approach, in his study of the 
asymptotic distributions of the poles of the integrale tritronquee. 

A priori, the WKB method is expected to give an approximation of poles 
z = a for a sufficiently large. Surprisingly our approach proves to be nu- 
merically very efficient also for poles close to the origin, see Table 2 below. 


The paper is organized as follows. In Section 2 we derive the Schrodinger 
equation associate with P-I and study thoroughly its relations with poles of 
P-I transcendents. Section 3 is devoted to the topological classification of 
Stokes complexes. In Section 4 we calculate the monodromy data in the 
WKB approximation, we derive the correct Bohr-Sommerfeld conditions for 
the poles of tritronquee, and we introduce the "small parameter" of the 
approximation. In Section 5 we obtain an asymptotic description of poles 
of the integrale tritronquee. In Appendix A and Appendix B we prove some 
theorems regarding the WKB functions that are used in section 2 and 3. 
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the theory of nonhnear waves and their appHcations". 



2 Poles and Cubic Oscillators 

We review some well-known facts regarding the isomonodromic approach to 
the P-I equation and analyze the isomonodromic deformation in a neighbor- 
hood of the singularities. 

2.1 P-I as an Isomonodromic Deformation 

P-I is equivalent to the compatibility condition of the following system of 
linear ODEs: 



-$jXz) - ( y'^'^ 2X^ + 2\y{z)-z + 2y\z)\^ 

^.(A,.) = ^^(^) + ^)t(A,.) . (2) 

The precise meaning of the word compatibility is given by the following 

Lemma 1. Fix zq, Aq and the Cauchy data y{zo), y'{zo), and <I>(Ao,^;o)- 
Let be any simply connected neighborhood of zq. Then y{z) satisfies the 
Painleve first equation in Uzq iff the system has a solution V(A, z) G 

C X Uzo- Moreover the solution is unique. 

Proof. See |Kap04| . □ 

In this subsection we suppose that we have fixed a solution y of P-I and 
a simply connected region U such that y\ij is holomorphic. 

We are now going to define the important concepts of monodromy data 
and isomonodromic deformation of equation ([1]). For this reason, we have to 
introduce some particular solutions of system (|l|2p . to be uniquely defined 
by the asymptotic behaviour for A — t- cxo. 

Fix A: G Z5 = {—2, . . . , 2} and the branch of A2 in such a way that 
ReAa — )• -|-oo as |A| — )• oo,argA = Then (see |Kap04| ) for any y 

solution of P-I, there exists a unique solution $fc(A,2;) of (|l|2p such that 




hm e+t-^-t-M- ■ Zi]^kiX)={]],yzeU, (3) 



A— ^00 



where A 4 is defined globally on the complex plane minus the negative real 
axis, and is positive on the positive real axis. Notice that, depending on 
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A: E Zs, (^A^j may not be equal to A2. Here and in the following, if not 
otherwise stated, e is an arbitrarily small positive number. 

From the asymptotics ([3| it follows that ^fc(A,z) and <1> fc_|_i(A, z) are 
linearly independent for any A; E Z5 and the following equality holds true 

^fc_i(A) = ^fc+i(A) + afc(z)^fc(A) , (4) 

where a^^z) is an analytic function of z, for any A; E Z5. 

Definition 1. Fixed z, we call the k-th Stokes multiplier of equation (QP 
and the set of all five Stokes multipliers the monodromy data of (Qp. The 
problem of calculating the monodromy data is called the direct monodromy 
problem. 

Stokes multipliers are very important for our analysis and we list their 
main properties in the following 

Lemma 2. Let ak{z),k E Z5 be defined as above. Then 

(i) equation ^ is an isomonodromic deformation of equation ([ip, i.e. 

dz ~ ^■ 

(a) The numbers ak,k G satisfy the following system of algebraic equa- 
tions 

1 + CTkak+i = -i crk+3 , A; E Z5 . (5) 
Proof See |Kap04| . 

□ 

Observe that only 3 of the algebraic equations (|5]) are independent. 

Definition 2. We denote V the algebraic variety of quintuplets of complex 
numbers satisfying ^ and call admissible monodromy data the elements of 
V. Due to Lemma\^ equations define the following map 

M : {P'l transcendents} — )■ V . 

Lemma. Ai is injective. 

Proof See |Kap04| . □ 

We end the section with a result of Kapaev, which completely character- 
izes the integrale tritronquee in term of Stokes multipliers. 

Theorem 1. (Kapaev) 

The image under M of the integrale tritronquee are the monodromy data 
uniquely characterized by the following equalities 

0"2 = cr-2 = . (6) 
Proof See |Kap04| . □ 
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2.2 Poles of y: cubic oscillator 

So far we dealt with the system (|l|2p in a region U which does not contain 
any pole of y{z). Indeed, the situation at a pole is different, for equation ([!]) 
makes no sense. _^ 

However, we show that any solution <I> (A, z) of system dH |2|) is mero- 
morphic in all the z — plane] moreover, a pole 2; = a of y{z) is also a pole 
of ^(A, z) and the residue at the pole of its second component satisfies the 
scalar equation of Schrodinger type (|8]). 

In order to be able to describe the local behavior of ^(A , z) near a pole 
a of we have to know the local behavior of y{z) close to the same point 
a. 

Lemma 3 (Painleve). Let a £ C be a pole of y. Then in a neighborhood of 
a, y has the following convergent Laurent expansion 



1 a{z - a)2 {z — a)^ 



+ ^^^±-^+'-^-^ + b{z-a)' + Y^c,{aMz-ay (7) 



(z - a)2 10 6 

j>5 



where b is some complex number and Cj{a,b) are real polynomials in a and 
b, not depending on the particular solution y. 

Conversely, fixed arbitrary a,b £ C, the above expansion has a non zero 
radius of convergence and solves P-L 

Proof See |GLSOO) . □ 

Definition 3. We define the map 

£ : — > {P-I transcendents} . 

C{a,b) is the unique analytic continuation of the Laurent expansion 

We have already collected all elements necessary to formulate the impor- 
tant 

Theorem 2. Fix a solution y of P-I and let <I>1*^(A, 2:), i = 1,2 /c S Z5 be 



the i-th component o/'^fc(A,z). Th 



k 
en 



(i) <I>fc(A, 2) is a meromorphic function of z. All the singularities are 
double poles. Moreover, a £ C is a pole of ^fc(A, z) iff it is a pole of y. 

(a) If a £ C is a pole of y then 

^k{X) = lim (z-a) $l'^(A,z) 

is an entire function of X. It satisfies the following Schrodinger equation 
with cubic potential 

^^^^ = (4A3 - 2aA - 286) ^k{X) , (8) 
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where b £ C is the coefficient entering into the Laurent expansion ^ 
of y around a. 

(Hi) If and A* are chosen as in asymptotics then Ve > 

lim A3e+t^^-^'^^^*fc(A) = i (9) 

' I 5~l ~5 

(iv) Equation ^ possesses the same monodromy data as equation i.e. 

^'fe_i(A) = ^'fe+i(A) + (Tfc^fc(A) . 

Proof, (i) From the Laurent expansion ([7|, it is easily seen that a pole a 
of y is a fuchsian singularity with trivial monodromy of equation ([2]). In 
particular the following Laurent expansions of ^f{\z) are valid 

^k\^^^) = ^^^{l-^{z-af^+^k{X){z-af+0{{z-af), 

^k\^^^) = ^^^(l + ^(^-a)2)-2(^,(A)(z-a) + O((z-a)2)(10) 

Expansions (jlOp show that $/fc(A,z) is meromorphic in a neighborhood of 
the point a and this point is a pole of order not greater than 2. 
(ii), (iii) The proof is in Appendix [Bl 

(iv) Since the functions {z — satisfy equations (|4]) for any z with 

constant Stokes multipliers, then their limit, i.e. the functions ^'^(A), satisfy 
the same equations. 

□ 

Definition 4. We call any cubic polynomial of the form V{X; a, b) = 4A^ — 
2a A — 286 a cubic potential. The above formula identifies the space of cubic 
potentials with 3 (a, b). 
We define the map 

r-.c'^^v. 

T{a,b) is the monodromy data of equation 

Theorem [2] has the following 
Corollary. M q C=T. 

The above corollary implies 

Theorem 3. Let y be any solution of P-L Then a £ C is a pole of y iff there 
exists 6 G C such that M{y) = T{cl, b). 

We finish this section with a theorem from Nevanlinna's theory |Nev70| . 
which implies the surjectivity of the map Ai. 



7 



Theorem 4. The map T is surjective. The preimage of any admissible 
monodromy data is a countable infinite subset of the space of cubic potentials. 



Proof See plf34] . □ 

As a consequence of the above theorem we have 

Theorem 5 (stated in |KK93| ). The map M is bijective: solutions of P-I 
are in 1-to-l correspondence with admissible monodromy data. 

Theorem [3] shows that the distribution of poles of P-I transcendents is 
a part of the theory of anharmonic oscillators^ which has been object of 
intense study since the seminal papers |BW68| and |Sim70| . 

Remark. In the theory of anharmonic oscillators a special importance is 
given to the vanishing of some Stokes multipliers. For a given k £ Z^, the 
problem is to find all (a, b) G such that the Stokes multiplier of equation 
^ vanishes. This is called the k-th lateral connection problem. Since fixed a, 
there exists a discrete number of solutions to any lateral connection problem, 
equation ak = is referred to as a quantization condition. 

As a consequence of Theorem [1] and Theorem |3l we have the following 

Corollary. The point a £ C is a pole of the integrale tritronquee if and 
only if there exists b £ C such that the Schrddinger equation with the cubic 
potential V{X; a, b) admits the simultaneous solution of two different quanti- 
zation conditions, namely a±2 = 0. 



2.3 Asymptotic values 

As it was previously observed, Stokes multipliers are defined by particular 
normalized solutions of equations ([TJ and ([8|. Following Nevanlinna, we 
define the monodromy data of equation ([8]) in a more invariant way. 

Definition 5. Let {ip, x} be a basis of solution of 
We call 

Wk{^,x)= lim ^ eCUoo, /cGZs . (11) 

the k-th asymptotic value. 

We collect the main properties of the asymptotic values in the following 

Lemma 4. (i) Let (p' = acp -\-bx o-nd x' = c(p -\- dx' , a,b,c,d £ C Then 

, , awkip,x) + b , . 

Wkif ,X) = 7 ^—7 • 12 
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(in) Wk+i{(p,x) / Wk{ip,x) 



Proof. See pf34] . 



□ 



Making use of equation (|lip . given the Stokes multipliers it is possible 
to calculate the asymptotic values. The converse is also true. In particular, 
the asymptotic values associated to the tritronquee integrale can be chosen 
to be 



3 Stokes Complexes 

In the complex WKB method a prominent role is played by the Stokes and 
anti-Stokes lines, and in particular by the topology of the Stokes complex, 
which is the union of the Stokes lines. 

The main result of this section is the Classification Theorem, where we 
show that the topological classification of Stokes complexes divides the space 
of cubic potentials into seven disjoint subsets. 

Even though Stokes and anti-Stokes lines are well-known objects, there 
is no standard convention about their definitions, so that some authors call 
Stokes lines what others call anti-Stokes lines. We follow here the notation 
of Fedoryuk |Fed93) . 

Remark. To simplify the notation and avoid repetitions, we study the 
Stokes lines only. Every single statement in the following section remains 
true if the word Stokes is replaced with the word anti-Stokes, provided in 
equation ^14\ ) the angles (pk are replaced with the angles ipk-\-^. 

Definition 6. A simple (resp. double, resp. triple) zero Xi of V{X) = 
V{X;a,b) is called a simple (resp. double, resp. triple) turning point. All 
other points are called generic. 

Fix a generic point Aq and a choice of the sign of ^/V{Xo). We call action 
the analytic function 



defined on the universal covering of X-plane minus the turning points. 

Let be the level curve of the real part of the action passing through 
a lift of Aq- Call its projection to the punctured plane ix^. Since i\^ is a 
one dimensional manifold, it is diffeomorphic to a circle or to a line. If i^o 
is diffeomorphic to the real line, we choose one diffeomorphism i\^^{x),x G M 
in such a way that the continuation along the curve of the imaginary part of 
the action is a monotone increasing function of x G M. 



Wq = Q,Wi = W-2 = 1,W2 = W-i = OO . 



(13) 
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Lemma 5. Let Aq be a generic point. Then i^o diffeomorphic to the real 
line, the limit lima;_)>+oo ^Aq (^) exists (as a point in C IJ ooj and it satisfies 
the following dichotomy: 

(i) Either liuix-^+ooixoix) = oo and the curve is asymptotic to one of the 
following rays of the complex plane 

A = pe^^^(^fc = ^^^^^,pGM+,A:GZ5, (14) 
5 

(a) or lim2:^.+oo ^Ao (^) = ^ij where Aj is a turning point. 
Furthermore, 

(Hi) i/lima;_j.±oo ^Ao(^) = oo then the asymptotic ray in the positive direction 
is different from the asymptotic ray in the negative direction. 

(iv) Let <pk,k E Z5 be defined as in equation ( |i^[ ). Then Ve > 0,3K E 
such that if (pk-i + e < argAo < <pk — £ and |Ao| > I^, then 
lima;_j.±oo ^Ao (^) = CO. Moreover the asymptotic rays of ix^ are the 
ones with arguments ipk and ipk-i- 

Proof. See |Str84]. 

□ 

Definition 7. We call Stokes line the trajectory of any curve ix^ such that 
there exists at least one turning point belonging to its boundary. 

We call a Stokes line internal if 00 does not belong to its boundary. 

We call Stokes complex the union of all the Stokes lines together with the 
turning points. 

We state all important properties of the Stokes lines in the following 

Theorem 6. The following statements hold true 

(i) The Stokes complex is simply connected. In particular, the boundary 
of any internal Stokes line is the union of two different turning points. 

(a) Any simple (resp. double, resp. triple) turning point belongs to the 
boundary of 3 (resp. 4, fesp 5) Stokes lines. 

(Hi) If a turning point belongs to the boundary of two different non-internal 
Stokes lines then these lines have different asymptotic rays. 

(iv) For any ray with the argument ipk as in equation ^14^ , there exists a 
Stokes line asymptotic to it. 

Proof See (StrM]. □ 
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3.1 Topology of Stokes complexes 

In what follows, we give a complete classification of the Stokes complexes, 
with respect to the orientation preserving homeomorphisms of the plane. 
We define the map L from the A-plane to the interior of the unit disc as 

L : Di 

Lipe'^P) = -e^^arctanp. (15) 
vr 

The image under the map L of the Stokes complex is naturally a deco- 
rated graph embedded in the closed unit disc. The vertices are the images of 
the turning points and the five points on the boundary of the unit disc with 
arguments (pk, with ipk as in equation (|14p . The bonds are obviously the 
images of the Stokes lines. We call the first set of vertices internal and the 
second set of vertices external. External vertices are decorated with the num- 
bers k £ Z5. We denote S the decorated embedded graph just described. 
Notice that due to Theorem [H] (iii), there exists not more than one bond 
connecting two vertices. 

The combinatorial properties of S are described in the following 

Lemma 6. S possesses the following properties 

(i) the sub-graph spanned by the internal vertices has no cycles. 

(a) Any simple (resp. double, resp. triple) turning point has valency 3 
(resp. 4, resp. 5). 

(iii) The valency of any external vertex is at least one. 

Proof. (i) Theorem |6] part (i) 

(ii) Theorem [6] part (ii) 
(iii) Theorem [6] part (iv) 

□ 

Definition 8. We call an admissible graph any decorated simple graph em- 
bedded in the closure of the unit disc, with three internal vertices and five 
decorated external vertices, such that (i) the cyclic-order inherited from the 
decoration coincides with the one inherited from the counter-clockwise orien- 
tation of the boundary, and (ii) it satisfies all the properties of Lemma\^ We 
call two admissible graphs equivalent if there exists an orientation-preserving 
homeomorphism of the disk mapping one graph into the other. 

Theorem 7. Classification Theorem 

All equivalence classes of admissible maps are, modulo a shift k ^ k -\- 
m,m £ of the decoration, the ones depicted in Figure [IJ 
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Figure 1: All the equivalence classes of admissible graphs. 
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Proof. Let us start analyzing the admissible graphs with three internal ver- 
tices and no internal edges. 

Any internal vertex is adjacent to a triplet of external vertices. Due to 
the Jordan curve theorem, there exists an internal vertex, say Aq, adjacent 
to a triplet of non consecutive external vertices. Performing a shift, they 
can be chosen to be the ones labelled by 0, 2, —1. Call the respective edges 
eo,e_i,e2. 

The disk is cut in three disjoint domains by those three edges. No internal 
vertices can belong to the domain cut by eo and 64, since it could be adjacent 
only to two external vertices, namely the ones labelled with and —1. By 
similar reasoning it is easy to show that one and only one vertex belong to 
each remaining domains. 

Such embedded graph is equivalent to the graph (300). 

Classifications for all other cases may proved by similar methods. 

□ 

The equivalence classes are encoded by a triplet of numbers (a b c): a is 
the number of simple turning points, b is the number of internal Stokes lines, 
while c is a progressive number, distinguishing non-equivalent graphs with 
same a and h. Some additional information shown Figure [1] will be explained 
in the next section. 

Remark. For any admissible graph there exists a real polynomial with an 
equivalent Stokes complex. 

Remark. Notice that the automorphism group of every graph in Figure [i] 
is trivial. Therefore the unlabelled vertices can be labelled. In the following 
we will label the turning points as in figure [TJ We denote "Boutroux graph " 
the graph (320) 

3.2 Stokes Sectors 

In the A-plane the complement of the Stokes complex is the disjoint union of 
a finite number of connected and simply-connected domains, each of them 
called a sector. 

Combining Theorem |6] and the Classification Theorem we obtain the 
following 

Lemma 7. All the curves ix^, with Aq belonging to a given sector, have the 
same two asymptotic rays. Moreover, two different sectors have different 
pairs of asymptotic rays. 

For any A; G Z5 there is a sector, called the k-th Stokes sectors, whose 
asymptotic rays have arguments (pk-i and ipk. This sector will be denoted 
Sfc. The boundary dT^k of each is connected. 
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Any other sector has asymptotic rays with arguments <pk-i and (fk+i, for 
some k. We call such a sector the k-th sector of hand type, and we denote it 
Bk- The boundary dBk of each Bk has two connected components. 

Choose a sector and a point Aq belonging to it. The function S'(Ao, A) is 
easily seen to be bi-holomorphic into the image of this sector. In particular, 
with one choice of the sign of ^/V it maps a Stokes sector into the half plane 
ReS > c, for some — oo < c < while it maps a B}^ sector in the vertical 
strip c < Re/S < d, for some — oo < c < < d < +00. 

Definition 9. We call a differentiahle curve 7 : [0, 1] — )• C an admissible 
path provided 7 is injective on [0, 1[, Aj ^ 7([0, 1]), for all turning points Aj, 
and ReS'(7(0), 7(t)) is a monotone function of t £ [0, 1]. 

We say that T,j ^ if there exist fij S Sj, fi^ G '^k and an admissible 
path such that 7(0) = /_fj,7(l) = fi^- 

The relation ^ is obviously reflexive and symmetric but it is not in 
general transitive. 

Notice that T,j ^ if and only if for every point fij £ Sj and every 
point fik G an admissible path exists. 

Lemma 8. The relation ^ depends only on the equivalence class of the 
Stokes complex S. 

Proof. Consider an admissible path from Sj to T,^, j 7^ k. The path is 
naturally associated to the sequence of Stokes lines that it crosses. We denote 
the sequence Imn = 0, . . . , A^, for some N € N. We continue analytically 
S{iJ,j, •) to a covering of the union of the Stokes sectors crossed by the path 
together with the Stokes lines belonging to the sequence. Since S{fij, •) is 
constant along each connected component of the boundary of every lift of a 
sector crossed by the path, then each of such connected components cannot 
be crossed twice by the path. Hence, due to the classification theorem no 
admissible path is a loop. Therefore, the union of the Stokes sectors crossed 
by the path together with the Stokes lines belonging to the sequence is simply 
connected. 

Conversely, given any injective sequence of Stokes lines /„, n = . . . ,N 
such that for any 0<n<A^— 1, /„ and In+i belong to two different 
connected components of the boundary of a same sector, there exists an 
admissible path with that associated sequence. This last observation implies 
that the relation ^ depends only on the topology of the graph S. Moreover, 
if the sequence exists it is unique; indeed, if there existed two admissible 
paths, joining the same fij and /ifc but with different sequences, then there 
would be an admissible loop. □ 

With the help of Lemma [8] and of the Classification Theorem, relation 
^ can be easily computed, as it is shown in Table [1] As it is evident from 
Figure dl for any graph type we have that ^ ^k+i, V/c S Z5. 
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Map 


Pairs of non consecutive Sectors not satisfying the relation ^ 


300 


None 


310 


(So, S2), (So, S_2) 


311 


(Si,S_i) 


320 


(Si,S„i), (Si,S_2), (S_i,S2) 


100 


(Si,S_i), (So,Il_2), (^0,^2) 


110 


All but (Si,S_i) 


000 


All 



Table 1: Computation of the relation ^ 



4 Complex WKB Method and Asymptotic Values 

In this section we introduce the WKB functions jk,k G Z5 and use them to 
evaluate the asymptotic values of equation (|8|). The topology of the Stokes 
complex will show all its importance in these computations. 
On any Stokes sector S^, we define the functions 



5fc(A) = 5(A*,A), (16) 
1 f'^ V'(y) 

j^(A) = e-^^(^)+-^'=(^) . (18) 

Here A* is an arbitrary point belonging to and the branch of VV is 
such that ReS'fc(A) is bounded from below. 
We call jfc the k-th WKB function. 

4.1 Maximal Domains 

In this subsection we construct the k-th maximal domain, that we denote 
-Dfc. This is the domain of the complex plane where the k-th WKB function 
approximates a solution of equation (|8|). 

The construction is done for any fc in a few steps (see Figure |2] for the 
example of the Stokes complex of type (300)): 

(i) for every S; such that ^ S^, denote D^^i the union of the sectors 
and of the Stokes lines crossed by any admissible path connecting 
and Sfc. 

(ii) Let Dk = Uz-^/,fc- Hence is a connected and simply connected 
subset of the complex plane whose boundary dD^ is the union of some 
Stokes lines. 
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(iii) Remove a (5-tubular neighborhood of the boundary dD^, for an arbi- 
trarily small 6 > 0, such that the resulting domain is still connected. 

(iii) For all I ^ k,l ^ k — 1, remove from Dk an angle A = pe*"^, \(f — (fi\ < 
e, p > R, for e arbitrarily small and R arbitrarily big, in such a way 
that the resulting domain is still connected. The remaining domain is 



4.2 Main Theorem of WKB Approximation 

We can now state the main theorem of the WKB approximation. 

Theorem 8 (G.D. Birkhoff |Bir33) . Olver |01v74) ). Continue the WKB func- 
tion jk to Dk- Then there exists a solution ipk{^) of such that for all 
XeDk 

< 5(A) (e2^(^) - l) 

)^(A)(e2pW_i) 

Here pk is a bounded positive continuous function, called the error func- 
tion, satisfying 

lim /5^(A)=0, 

A— >oo 
V3fe_i<argA<95fe_,_i 

and g{X) is a positive function such that g{X) < 1 and 

lim g{X) = J . 

A— >oo Z 

AeDfenEfe±2 

Proof. The proof is in the appendix [A] □ 

Notice that is sub-dominant (i.e. it decays exponentially) in and 
dominant (i.e. it grows exponentially) in Si,V/ 7^ k. 

For the properties of the error function, ■0^ is subdominant in and 
dominant in T,k±i- Therefore, in any Stokes sector there exists a sub- 
dominant solution, which is defined uniquely up to multiplication by a non 
zero constant. 



V^fc(A) 
jk{X) 



+ 1 



< 



V'iX) 



Av{xy. 



+ (1 + 



V'iX) 



Av{xy. 
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Figure 2: In the drawings, the construction of Dq for a graph of type (300) 
is depicted. 



17 



4.3 Computations of Asymptotic Values in WKB Approxi- 
mation 

The aim of this paragraph is to compute the asymptotic values for the 
Schrodinger equation (j8]) in WKB approximation. We explicitely work out 
the example of the Stokes complex of type (320), relevant to the study of 
poles of the integrale tritronquee. 

Definition 10. Define the relative errors 



lim pi{X), if ^ Sfc 

A— >oo 
AGSfcHA 

oo, otherwise 



and the asymptotic values 



Wk{l,m)'^= Wkiipui^m)- (19) 



We say that ~ provided pf < The relation ^ is a sub-relation of 

Notice that p^^^ = and = p^^ (see Appendix IX]) . 
In order to compute the asymptotic value wtil, m), we have to know the 
asymptotic behavior of ipi and in S^. By Theorem [8l 



A^-oo ?';(A) 



lim 



/O, if ^(e^/'i - 1) < 1 . 



Hence the asymptotic behavior of ipi in can be related to the asymptotic 
behavior of ji in if the relative error pf is so small that the above inequality 
holds true, i.e. if S/; ~ E;. 

Remark. Depending on the type of the graph S, there may not exist two 
indices k ^ I such that all the relative errors pf,p'^,n £ Z5 are small. How- 
ever it is often possible to compute an approximation of all the asymptotic 
values Wn{l,k) using the strategy below. 

(i) We select a pair of non consecutive Stokes sectors S;, S;_|_2, with the 
hypothesis that the functions -0; and '0i+2 ^-re linearly independent, so 
that wi{l, 1 + 2) =0, wi+2{l, 1 + 2) =00. Since p\~^^ = p\'ll = then 

Wi+i{l,l + 2)= lim ^^^^^ 



A-s>oo 7';+2(A) 

AeE,+inDinA+2 

Therefore, we find three exact and distinct asymptotic values. 



18 



(ii) For any k ^ 1,1 + 1,1 + 2 such that ~ and ~ S^, we define 
the approximate asymptotic value 



Wk{l,m) 



lim 



jm{z) 



The spherical distance between Wk{l,m) and wi.(l,m) may be easily 
estimated from above knowing the relative errors and p^^^- 

If for any k ^ 1, 1 + 1,1 + 2, ~ and ~ S^, then the calculation 
is completed. 

(iii) If not, we can use the fact that quintuplets of asymptotic values for 
different choices of ipi , ipi+2 are related by a Mobius transformation 
(see formula (|12p ). If for some pair + 2) the assumption ~ 
Sfc, ~ Sfc fails to be true for just one value of the index k = k* , 
and, for another pair (/',/' + 2) the assumption T^i ~ S^/, ~ 
Sfc/ fails to be true for just one valued of the index k' = k'*, with 
k'* 7^ k* , then there are three values of the index m € Z5 such that an 
approximation of Wm{l,l + 2) and Wmil',l' + 2) is computable. Since 
any Mobius transformation is fixed by the action on three values, then 
we can compute an approximation of the transformation relating the 
quintuplets Wk{l,l + 2) and Wk{l' , I' + 2) for any S Z5. Hence, we can 
calculate an approximation of the whole quintuplets Wk{l,l + 2) and 



Remark. As shown in TableUl the relation ^ is uniquely characterized by 
the graph type. For the sake of computing the asymptotic values the important 
relation is ~ and not Indeed, the calculations for a given graph type, say 
(a b c), are valid for (and only for) all the potentials whose relation ~ is 
equivalent to the relation ^ characterizing the graph type (aba). 

Due to the above remark, in what follows we suppose that the relation 
~ is equivalent to the relation We have the following 

Lemma 9. Let V{\;a,b) such that the type of the Stokes complex is (300), 
(310), (311); moreover, suppose that the ~ relation coincides with Then 
all the asymptotic values of equation ^ are pairwise distinct, but for at most 
pair. 

Proof. For a graph of type (300) or (311) the thesis is trivial. For a graph 
of type (320), it may be that wq = W2 or wq = W-2. Since W2 7^ W-2 the 
thesis follows. □ 

We completely work out the case of Stokes complex of type (320), while 
for the other cases we present the results only. Due to Lemma |9l we omit 
the results for potentials whose graph type is (300), (310) and (311). 



:',/' + 2). 
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Boutroux Graph = 320 We suppose that Sq ~ '^±2- 

Let us consider first the pah Sq and S_-2- In Figure [3] the maximal do- 
mains Dq ans D-2 are depicted by colouring the Stokes lines not belonging to 
them blue and red respecitvely. In particular SQ,Lo,jQ (resp. S-2,L-2,j-2) 
can be extended to all Dq (resp. -D_2) along any curve that does not intersect 
any blue (resp. red) Stokes line. 




Figure 3: Calculation of W-i{0, —2) and of W2{0, —2) 

We fix a point A* G Sq such that So{X*) = S'_2(A*) = Lo(A* 
L_2(A*) = 0. 
By definition 

w,{0,-2) = lim 



A->oofc j-2(A) 

= lim e-'^o(^)+'^-2(^)e^"(^)--^-2(-^) 

X—^oOf^ 

Here A — t- ook is a short-hand notation for A — t- oo, A € Syt n -Dq ^ D^2- We 
calculate Wk^O, —2) for k = —1, 2. 

We first calculate limA^oo^ e-'^o(^)+^-2(A)^ 

Notice that ^ = in S^.. Hence 

lim -5o(A) + 5_2(A) = -Soink) + 5'_2(/ifc) , A; = -1, 2 , 

A— >oOfe 

where //^ is any point belonging to (in Figure |3l the paths of integration 
defining 5'o(/ifc) and S-2{fik) are coloured blue and red respectively). 
On the other hand, since = — ^q^^ in Sq |J 5]_2, we have that 

—SoifJ^k) + S-2{fJ'k) = — 25o(As) , s = — 1 if A; = — 1 and s = if fc = 2 . 
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We now compute limA^oo,. e^^^^^-^-^^^). Since = in f] L>_2, 
we have that 

hm Lo(A) - L_2(A) = -Lo(Atfe) - i-2(/Ufc) , 

A— >oOfc 

Lo(^fc)-L-2(R0 = ^^4^,/c = -l,2. 

Here is the blue path connecting A* with composed with the inverse 
of the red path connecting A* with /i^ (see Figure |3]) . 
Therefore, we have 

2-Ki 

hm -Lo(A) — L_-2(A) = — — , o" = —1 if /c = —1 and o" = +1 if A; = 2 . 

A— >oOfc 4 

Combining the above computations, we get 

We stress that W-i{Q, —2) is exact while W2{Q-, —2) is an approximation. 
Performing the same computations for the pair Sq and S2, we obtain 



wi{<d, 2) = -i e-2^o(^i), -w.aCO, 2)=ie 



■25o(Ao) 



Having calculated the triplet of asymptotic values wo,W2,W-2 for two 
different pairs of Stokes sectors, we can compute an approximation of the 
Mobius transformation relating all the asymptotic values for the two pairs: 

We eventually compute the last asymptotic value for the pair So,S_2, 
that is 



p-25o(Ai) 

w^i(0,-2) = 



1 + e-2('5o(Ai)-So(Ao)) 



Quantization Conditions The computations above provides us with 
the following quantization conditions: 

= w^2 ^ e~2('^o(^i)~'^o(^o)) = -1 (20) 

= ^i;„l4^e-2(^o(^-)-^o(^«)) = -l (21) 
= ^_^^e-2(So{Ai)-So(A_i)) = _l + e-2(So(Ai)-5o(A„)) (22) 



We notice that equation (|22p is incompatible both with (|20p and (|2ip . 
Equations (|20p and (j2T|) are Bohr-Sommerfeld quantizations. 
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As was shown in equation (|13p . the poles of the integrale tritronquee 
are related to the polynomials such that wi = w^2 and w^i = W2- Since 
equations (|20p and (|2ip can be simultaneously solved, solutions of system 
(|20|2ip describe, in WKB approximation, polynomials related to the in- 
tegrale tritronquee. System (|20|2ip was found by Boutroux in |Boul3| 
(through a completely different analysis), to characterize the asymptotic dis- 
tribution of the poles of the integrale tritronquee. Therefore we call (|20|2ip 
the Bohr-Sommerfeld'Boutroux system. 

Equation (|22p will not be studied in this paper, even though is quite 
remarkable. Indeed, it describes the breaking of the PT symmetry (see 
[DTOO] and |BBM+01) V 

Case (100) 

«;o(l,-l) = -1 
uJ_2(l,-l) = W2{l,-l) = l 

Since wq 7^ u]±2 and W2 7^ W-2-, if the error or p^_^ is small enough, 
then all the asymptotic values are pair wise distinct. 

Case (110) 

«}_i(l,-2) = 1 
W2{l,-2) = -1 

In this case, it is impossible to calculate wq with the WKB method that 
has been here developed. Hence it may be that either wq = W2 or u^o = W-2- 

Notice, however, that (110) is the graph only of a very restricted class of 
potentials namely V{\) = (A -|- Ao)^(A — 2Ao), where Aq is real and positive. 
Since the potential is real then wq 7^ w±2- 

Case (000) In this case, no asymptotic values can be calculated. Notice, 
however, that V[X) = A^ is the only potential with graph (000). For this 
potential the asymptotic values can be computed exactly, simply using sym- 
metry considerations. Indeed one can choose Wk = e~* , k G Z5. 

4.4 Small Parameter 

The WKB method normally applies to problem with an external small pa- 
rameter^ usually denoted h or e. In the study of the distributions of poles of 
a given solution y of P-I there is no external small parameter and we have 
to explore the whole space of cubic potentials. The aim of this section is to 
introduce an internal small parameter in the space of cubic potentials, that 
greatly simplifies our study. 
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On the linear space of cubic potentials in canonical form 

V{X;a,b) = 4A^ - 2aX - 286, 

we define the following action of the group M"*" x Z5 (similar to what is 
called Symanzik rescaling in |Sim70| ) 

{x,m)[V{X;a,b)] = V{\;n^"'x'^a,n^'^x^b), x G M+, m G Z5, = (23) 

The induced action on the graph S, on the relative error pY^, and on the 
difference Si{Xj) — Si{Xk) is described in the following 

Lemma 10. Let the action of the group x Z5 be defined as above. Then 

(i) {x,m) leaves the graph S invariant, but for a shift of the labels k — >• 
k + m of the external vertices. 

(a) {x,m)[Si{Xj) - Si{Xk)] = X2 {Si{Xj) - Si{Xk)) ■ 

(Hi) {x,m)[p^] = x~2p'^ . 

Proof. The proof of (i) and (ii) follows from the following equality 

^ViX;n'^''x^a,n3''x3b)dX = xi y/ViX';a,b)dX' , X = xX! . 

The proof of point (iii) follows from a similar scaling law of the 1-form 
a{X)dX (see equation (|3ip in Appendix |A]) . □ 

Due to Lemma [TOT iii). e = ||| plays the role of the small parameter. 
Indeed, along any orbit of the action of the group M x Z5, all the (finite) 
relative errors go to zero uniformly as ||| — >• 0. 

Since all the relevant information is encoded in the quotient of the space 
of cubic potentials with respect to the group action, we define the following 
change of variable 

6 62 
u{a,b) = -, p{a,b) = ^ . (24) 

The induced action on these coordinates is simple, namely 

{x,m)[p{a,b)] = p{a,b) and (x, m)[zv(a, 6)] = Q"^x i'{a,b) . 

Moreover, the orbit of the set /i) € s.t. \u\ = 1, |argz^| < f , /x 7^ 0| 
is a dense open subset of the space of cubic potentials. 
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5 Poles of Integrale Tritronquee 



From Lemma |9] and the results of the computations in Section [4. 3| it follows 
that equation (|8|) admits in WKB approximation the simultaneous solutions 
of the two quantization conditions w±i = w^2 only if the Stokes complex is 
of type (320). In particular, after our calculations the poles of the integrale 
tritronquee are related, in WKB approximation, to the solutions of the Bohr- 
Sommerfeld-Boutroux system (|20l2ip . 

We rewrite this system in the following equivalent form: 



ai 



\/F(A; a, h)d\ = i7r(n — -) 



(25) 



^/V{\;a,h)d\ 



-i7r(m ) 

^ 2^ 



where m, n are positive natural numbers and the paths of integration are 
shown in figure H] 



branch cuts defining the 
■ square root of the potential 



(320) 




Figure 4: Riemann surface fj? = V{X; a, b) 

System (|25p is studied in detail in |KK93| where the following lemma is 
proven. 

Lemma 11. // a polynomial V{X) = 4A^ — 2aX — 28b satisfies the system 
(ED then |arga| > ^. 

The Lemma above should be compared to the conjecture, to which we 
referred in the introduction. 

Real Poles We compute all the real solutions of system (|25p and compare 
them with some numerical results from |JK01| . We note that the accuracy of 
the WKB method is astonishing also for small a and b (see Table 2 below). 
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For the subset of real potentials, we have 



J) yJV{\;a,b)d\ = I yJV{\]a,h)d\ , 

where ~ stands for complex conjugation. 

Therefore system (|25p reduces to one equation and the real poles of 
tritronquee are characterized, in WKB approximation, by one natural num- 
ber. 

Lemma 12. Let and v he defined as in equation ^24\ )- Then the real 
polynomials whose Stokes complex is of type (320) are the orbit of a single 
point of the action, characterized by fi* = —3158, 92 and v > 0. 

Moreover, if a real polynomial V{X; a, b) satisfies the Bohr-Sommerfeld- 
Boutroux conditions H25\) . then 

a = a*{n-l)l, b = b*{n-\)l (26) 
for some n G N* and a* ^ -4, 0874, b* ^ -0, 1470. 

Proof. All real cubic potentials whose Stokes complex has type (320) have 
one real turning point Aq and two complex conjugate turning points X±i. 
For the subset of real potentials with two complex conjugate turning points 
the cycles a±i, as shown in Figure |4l are unambiguously defined. From 
the Classification Theorem, it follows that a Stokes complex has type (320) 
if and only if Re£^ ^yV{X;a,b)dX = 0,lm£^ y/V{X; a,b)dX ^ 0. Since 
these conditions are invariant under the M"*" action, if they are satisfied for 
a point of the space of cubic potentials, then they are satisfied on all its 
orbit. Moreover, it is easily seen that this orbit exists and is unique. With 
the help of a software of numeric calculus, we characterized numerically the 
orbit. Afterthat, using the scaling law in Lemma [TOT ii) we calculated all the 
real solutions of the Bohr-Sommerfeld-Boutroux system. □ 

To the best of our knowledge, the asymptotic for the b coefficients has 
never been given. 

In the paper |JK01| , the authors showed that the integrale tritronquee has 
no poles on the real positive axis. The real poles are a decreasing sequence 
of negative numbers and some of them are evaluated numerically in the 
same paper. 

In Table 2, we compare the first two real solutions to system (|25p with 
the numerical evaluation of the first two poles of the integrale tritronquee. 

6 Concluding Remarks 

We have studied the distribution of the poles of solutions to the Painleve 
first equation using the theory of the cubic anharmonic oscillator. We have 
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WKB 


Numeric 


Error % 




-2,34 


-2,38 


1,5 


hi 


-0,064 


-0,062 


2 


/^i 


-3158 


-3510 


10 


a2 


-5,65 


-5,66 


0,2 


62 


-0,23 


unknown 


unknown 




-3158 


unknown 


unknown 



Table 2: Comparison between numerical and WKB evaluation of the first 
two real poles of the integrale tritronquee. 

applied a suitable version of the complex WKB method to analyze the dis- 
tribution of poles of the integrale tritronquee. 

In subsequent publications we plan to pursue our study of poles of P-I 
transcendents in different directions. 

In particular, we want to use the Nevanllina theory of the branched 
coverings of the sphere to complete the analysis of the poles of the integrale 
tritronquee, showing that the developed WKB method yields a complete 
qualitative picture and efficient quantitative estimates of the distribution of 
the poles. 

Since the quantization condition uo = characterizes the monodromy 
data of a family of special solutions of P-I, called integrale tronquee (see 
|Kap04| ), these solutions are strictly related to the spectral theory of FT sym- 
metric anharmonic oscillators and to functional equations of Bethe Ansatz 
type (for what concerns the PT symmetric anharmonic oscillators and the 
Bethe Ansatz equations, see |DDTOlj and references therein). We are going 
to investigate the consequences of this relation in a subsequent publication. 

Remark. After the main computations of the present paper had been com- 
pleted, the author learned from B. Dubrovin about the results of V. Novok- 
shenov presented at the conference NEEDS09 (May 2009). Novokshenov 
studied WKB solutions to the Schrddinger equation ^ with 6 = and their 
connections to the distributions of poles of certain particular solutions to the 
Painleve-I equation, including integrales tronquee and integrale tritronquee. 

Remark. In the paper JMBlOi/ written after this paper was published, the 
author proves that eventually (for big enough n and m ) around any solution 
of the Bohr-Sommerfeld-Boutroux system \25\) there is one and only one 
pole of the integrale tritronquee. Moreover the distance between a pole and 
its approximation vanishes asymptotically. 
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A Appendix 



The aim of this appendix is to prove Theorem [S] Our approach is similar to 
the approach of Fedoryuk |Fed93| . 

Notations are as in sections |3] and 21 except for oofc. In what follows, 
we suppose to have fixed a certain cubic potential ^(A; a, b) and a maximal 
domain D^- To simplify the notation we write V{X) instead of V{X; a, b). 

A.l Gauge Transform to an L-Diagonal System 

The strategy is to find a suitable gauge transform of equation ([8| such that 
for large A it simplifies. We rewrite the Schrodinger equation 

- <(A) + V{X)^P{X) = , (27) 

in first order form: 



*'(A) = ^(A)^'(A) 

' 1 
y(A) 



(T.?^^ n) • (28) 



Lemma 13 (Fedoryuk). In 

(i) the gauge transform 

Y(X) = A(X)U{X), 

AW = A(A)(^Lva« , (29) 

is non singular and 
(a) the system i28\) is transformed into the following one 

U'{X) = F{X)U{X) = {A-^EA - A~^A') U , 

Fix) = 2ynA)(j ;;)+o(.)(_\ _\), 



(30) 



«(A) = =5 (4m^"(A) - 5y'2(A)) . (31) 

32yy(A) 



Proof (i) Indeed det A{X) = 2jl{X)^/V{X) / 0, VA G Dk, by construc- 
tion of jyfc and D/j. 

(ii) It is proven by a simple calculation. 

□ 
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A. 2 Some Technical Lemmas 

Before we can begin the proof of Theorem [8l we have to introduce the com- 
pactification of and the preparatory Lemmas [14] and [T5j 

Compactification of Dk Since Dk is simply connected, it is conformally 
equivalent to the interior of the unit disk D. We denote U the uniformisation 
map, U : D ^ D^- 

By construction, the boundary of is the union of n free Jordan curves, 
all intersecting at oo. Here n is equal to the number of sectors S/ such that 
Efc minus 2. 

Due to an extension of Caratheodory's Theorem ( |Car| . §134-138), the 
map U extends to a continuous map from the closure of the unit circle to 
the closure of D/^. The map is injective on the closure of D minus the 
n counterimages of oo. Hence, the uniformisation map realizes a n point 
compactification of D^, that we call Df^. In there are n point at oo. We 
denote oo^ the point at oo belonging to the closure of U U U Tjk+i)- 

Moreover, for A = A: + 2 or A = A: — 2, if ^ we denote oo; the point at 
oo belonging to the closure of C/(S/). 

Definition 11. Let H be the space of function holomorphic in Dk and con- 
tinuous in Dk- H endowed with the sup norm is a Banach space (H , \\-\\h)- 
Let r(A),A G Dk — oofc he the set of injective piecewise differentiahle 
curves 7 : [0, 1] — t- Dk, such that 



4- the length of the curve restricted to [0,T] is O {\j{T)\) , as t — )• 1. 

Let r(A) be the subset ofT{X) of the paths along which ReSk{'y{0),'y{t)) 
is non decreasing. 

Let Ki : H ^ H and K2 : H ^ H be defined (for the moment formally) 



2. 



1 



7(0) = A, 7(1) = oofc, 

ReSk{'y{0),'y{t)) is eventually non decreasing, 



3. 



there is an e > such that eventually |arg7(i) — < ^ 




(33) 



(32) 



Let p: Dk^ Dk: 
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Remark. Since along rays of fixed argument if, with \(f — < ^ — 
ReSk is eventually increasing, there are paths satisfying point (1) through 
(4) of the above definition. Moreover, by construction of Dk, r(A) is non 
empty for any X. 

Before beginning the proof of the theorem, we need two preparatory 
lemmas. 

Lemma 14. Fix e > 0, an angle \aig ip — ^| < f ~ ^; let = n 
{A G C, |A - ^1 < f - e). Denote i{R) = i^^.^ nn,Re M+, and let L{R) 
be the length with respect to the euclidean metric ofi{R). Then L{R) = 0{R) 
and infygj(j^) |A'| = 0{R). 

Let r be any level curve of Si{X*, ■) asymptotic to the ray of argument 
n{R) = {A e n, ReSi{X, Re^'f') > O} , and M{R) be the length of r n n{R). 
Then M{R) = 0{R). 

Proof. |Str84| . chapter 3. □ 
Lemma 15. (i) p is a continuous function, 
(a) Ki and K2 are well-defined bounded operator. In particular 

mh]W\<p{X)\\h\\H,i = 1^2 (34) 

(ill) K2[h]{^k) = Ki[h]{^k) = Ki{<xk±2) = 0,yheH 

_ 7 

Proof. (i)Since a(A)ciA = 0(|A| '^), then a{X)dX is integrable along any 
curve 7 G r(A). Therefore p is a continuous function on D^- 

(ii)We first prove that (a) Ki[h]{X) does not depend on the integration 
path for any A € minus the points at infinity. A result that easily implies 
that Ki[h]{-) is an analytic function on Dk, continuous on A € Dk minus the 
points at 00. We then prove (b) the estimates (|34p and (c) the existence of 
the limits Ki[h]{ooi), I = oofc, 00^-1-2 • 

To simplify the notation, we prove the theorem for the operator Ki. The 
proof for K2 is almost identical . 

(a) Let ja^Jb S r(A). The curve i'y^(T): where T = 1 — e for some 
small e > intersect 7^ at some 7b(T'). Therefore we can decompose 
~lb ° 7a into two different paths with the help of a segment of ^-^^(t)) 
I-m e^^''^>'^^'>aip)h{p)dfi = J^^ + /^^ e2^'=(^'^)a(/i)h(/i)d/x. One path 71 
is the loop based at A and the other 72 is the loop based at ook- Since 
71 C Dk, then / ^ e^^^'^^'^^'i a{p)h{iJi)diJL = 0. Along 72, e2'5fc(72W,A) < ^ ^^lem- 
fore the integrand can be estimated just by |a(72(t))|. Due to lemma [T4| 

= 0(|7a(T)|~2 Since e is arbitrary, then Ki{h]{X) does not 
depend on the integration path. 

(b) Clearly for any path 7 E f(A), \Ki[h]{X)\ < /J \j{X)h{X)d)\ dt. Since 
Ki[h]{\) does not depend on 7, then estimate (j34|) follows. 



29 



(c) Let A„ be a sequence converging to ooi,l = k + 2 or I = k — 2; 
without losing any generality we suppose that the sequence is ordered such 
that ReSk{Xn) < -^^^^(An-i-i). Fix a curve r, as defined in Lemma [TU 
By construction of D^, it is always possible to connect two points A„ and 
Xn+m with a union of segments of the curves ix^,i\„^^ and of r. We de- 
note by 7 the union of this three segment. By construction of (see 
Subsection 14.11 (iii)), there exists e > such that |arg — ^| < f ~ e, Vn. 



+ |An+m|. Hence 
. Then K^[h]{Xn) 



Therefore, due to Lemma [141 7 has length of order |A. 

\K,[h]{Xn) - K,[h]{Xn+m)\ < \h{X)a{X)dX\ = 0{\Xr 
is a Cauchy sequence and the limit is well defined. 

We now prove that this limit is zero by calculating it along a fixed ray 
A = xe*"^ inside Sfc-|-2- Let us fix a point x* on this ray in such a way that the 



function ReSk{x* 
the ray we have 



is monotone decreasing in the interval [x* , +00 [. Along 



Ki[h]{x) 



ff e''^''(-y^'^'^a{y)h{y)dy + g{x* 



where g{x*) is a constant, namely /^gp(2,*) e'^^''^^''^* ^ a{fi)h{fi)dfj,. Hence 

lim^^oo Ki[h]{x) = lima;_^oo "'^f^'-''^ = 0. 

With similar methods the reader can prove that the limit Ki[h]{ook) 
exists and is zero. 



□ 



We are now ready to prove Theorem |5J 



Theorem 9. Extend the WKB function to D^. There exists a unique 
solution of ^ such that for all X G 



V'fc(A) 



ife(A)yy(A) 



1 



+ 1 



< g{X){e'PW-l) 



< 



V'{X) 



+ (1 + 



V'{X) 



MA)(e 



2p(A) 



where g{X) is a positive function, g{X) < 1 and g(oOfc-|-2) 



2- 



Proof. We seek a particular solution to the linear system (|30p via successive 
approximation. 

If U{X) = C/(^) ® U^'^'^ G if © if satisfies the following integral equation 
of Volterra type 



UiX) 
K[U]{X) 



Uo + K[U]{X) , Uo 



K2[C/« + [7(2)](A)y' 



(35) 
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then U{X) restricted to Dk satisfies (|30p . 
We define the the Neumann series as follows 



n+l 

Un+l = U'^ + K[Un],Un+l=Y,K'[U^]. (36) 

More explicitly, 



a(^2)(l - e2^(^2'^i)) . . . a(^„)(l - e^-^l^n.^n-i)) 
a(^2)(l - e2^(^2,Mi)) . . . a(^„)(l - e2^(^"'^"-i)) 

Here the integration path 7 belong to r(A). For any 7 G r(A) and any 
n > 1 



where i^"[C/o]'-*^ is the i-th component of ^"[C/o]- Hence 

\K-[UU\)\<~{2p{X)r (37) 

Thus the sequence ^7" converges in H and is a solution to (|35p ; call [/ 
its limit. Due to Lemma [T5l C/*'^-*(oofc±2) = 0. 

Let 'ffc be the solution to (j28)) whose gauge transform is [/ restricted to 
Dk] The first component ^/^^ of satisfies equation (|27p . 

From the gauge transform (|29p . we obtain 



^-1 = u,[X) + U2{\)-l, 



Jfc(A)VVXX) 4y(A)i 4y(A)t 

4y(A)i ' 

The thesis follows from these formulas, inequality (|37p and from the fact 
that f/i(oOfc_|_2) = 0. 

□ 
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Remark. The solution ipkW of equation ^ described in Theorem IB\ may 
be extended from to the whole complex plane, since the equation is linear 
with entire coefficients. The continuation is constructed in the following 
Corollary. 

Corollary 1. For any A G C, A not a turning point, we define T{\) as in 
Definition [771 Fixed any 7 G r(A) and h a continuous function on 7, we 
define the functionals Ki[K\{\) as in equations H32\) and i33\) . We define the 
Neumann series as in equations I\cl5\) and l\3b]) . and we continue jj. along 7. 

Then then Neumann series converges and we call U^'^\X) and U^'^\X) 
the first and second component of its limit. 

Moreover, V'fc(A) = (U^^\X) + C/^^^(A)) jfc(A) solves equation ^ and for 
any e > 

lim fc/W(A) + [/(2)(A)) = 1 

The reader should notice that if A ^ D^, then r(A) is empty and we 
cannot estimate ^.''i^) . 



B Appendix 

The aim of this Appendix is to prove Theorem |2] (ii) and (iii). The notation 
is, if not otherwise stated, as in the previous sections of the paper. 

Next to a pole z = a of a solution y{z) of P-I, equation ([1]) becomes 
meaningless. To get rid of this singularity we perform a gauge transform of 
([1]) such that the gauge-transformed equation has full meaning in the limit. 
In what follows, we suppose that z belongs to a punctured neighborhood of 
a, where y{z) is holomorphic. 

A gauge transform 

Let z be a fixed regular value of y{z). Let "^(A, z) = G(A, z)"^(A, z), 



I 2/'W+ 2(a4(.)) 



G(A,z) 



/ 



(38) 



\^2{\-y{z)) 

Then $ (A; z) satisfies ([1]) if and only if z) satisfies the following 

equation 



^(''^)=(q(A;z) 
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where 



Q{\;z)=A\'-2\z+2zv{z)-Ay%z)+y'\z)+ y'^''] + " (39) 

\-y{z) 4.[\-y{z)Y 

We denote '0 the first component of ^' . The equation for is equivalent 
to the following second order scalar equation for ip 



i'xx{X,z) = Q{X;z)iP{X,z) (40) 

We summarize some property of the perturbed potential, which can be 
easily verified using the expansion ([7|. 

Lemma 16. Let = = {z — a)^ + 0{{z — a)^) then 

(i) Q{\]z) has a double pole at X = It is an apparent fuchsian singu- 
larity for equation ( |^0[ ); the local monodromy around it is —1. 

(a) Q{X; z) has two simple zeros at X = + O(e^) 

(ill) Q{X;z) =4A3-2(a + e)A-286 + 0(e)-|^ + 4p^4^, where 0{e) 
does not depend on X. 

Equation (|40| is a perturbation of the cubic Schrodinger equation (|8]) 
and the asymptotic behaviours of solutions to the two equations are very 
similar. Indeed the local picture around the point at oo depends only on the 
terms 4A^ and —2zX^. 

More precisely, the equivalent of Corollary [1] in Appendix |A] is valid also 
for the perturbed Schrodinger equation. 

Definition 12. For any z, define a cut from X = (^^\y2 to oo such that it 
eventually does not belong to the the angular sector |argA — < 

Fix X* in the cut plane. Sk{X;z) = J^, \J Q{n; z)dfi is well-defined for 
|arg A — ^1^1 < ^ and X >> 0. Here the branch of \fQ is chosen such that 
ReSk{X) — ^ +00 as |A| CO, |argA - ^| < | - e. We define jk{X;z) as 
in equation H18\) and a{X; z) as in equation h31\) . but replacing V{X) with 
Q{X;z). 

For any X in the cut plane, let T{X) be the set of piecewise differentiable 
curves 7 : [0, 1] to the cut plane, 7(0) = A, 7(1) = 00, satisfying properties 
(2) (3) and (4) of Definition\M 

For any 7 G r(A), let H be the Banach space of continuous functions 
on 7 that have a finite limit as t ^ 1. Formulae l\32\) and l\33\) define two 
bounded functionals on H . We call such functionals Ki{X; z) and K2{X; z). 

Following the proof of Theorem [8l the reader can prove the following 
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Lemma 17. Let X belong to the cut plane, A not a zero of Q{-;z). Fixed 
any 7 G r(A), we define the Neumann series as in equations I\cl5\) and l\36]) . 
and we continue jk along 7. 

Then the Neumann series converges and V'fc(A) = (C/i(A) + C^2(A)) jfc(A) 
solves equation Moreover, for any e > 

lim fc/W(A) + [/(2)(A)) = 1 

Definition 13. Let i/^fe(A, z) be the unique solution of equation d^Op such 
that 

ibkiX, z) , , , , 2Trk 

^''^ ' ^ ' 1, as \X\^ 00, argX= — . (41) 



A-ie-s^^+^^^ 



3 4 , 



ifere i/ie branch of X^ is fixed as X ^ 00, argA = and there it 

coincides with the branch chosen in equation We define ^fc(A,a) to be 

the unique solution of equation ^ with asymptotic ( |^ j[ ), where z = a. 

We denote ipk{X;z) the unordered pair |'(^fc(A, z), — ^/c(A, 

Remark. We notice that if the cuts are continuous in z, then ipkiX,z) = 
c{z)'ipk{X), where ipkW is the solution constructed in Lemma\Tl\and c{z) is 
a bounded holomorphic function. 

Theorem 10. lim^_>,a VfcC-^) -2) = i^ki^^o), VA G C. 

Proof. Let A be any point in tlie complex plane which is not a zero of 
V{X;a,h). For any sequence e„ converging to zero, we choose two fixed 
rays ri and r2 of different argument ipi and ip2-, \pi — < f • We denote 
Dij^e a disk of radius R with center A = ^ and we split the sequence e„ into 
two subsequences such that rj Pi = for any n big enough. 

For i = 1,2, we choose the cuts defined in Definition [12] in such a way 
that there exists a differentiable curve 7j : [0, 1] — >■ C, 7i(0) = A, 7i(l) = 00 
with the following properties: (1)7^ avoids the zeroes of Q(A,e^) and a fixed, 
arbitrarily small, neighborhood of the zeroes of y(A;a,5), (ii)7j does not 
intersect any cut, and (iv) 7j eventually lies on r^. 

The proof of the thesis relies on the following estimates: 



SUp;,gC-DH,, 



A-^ 



A-^ 



A^ 



|Q(A; a + e) - V{X; o, 6)| = ©(e^'^-S) , (42) 
\Qx{X-,a + e)-Vx{X-,a,h)\=0{e^'-^) , 
\Qxx{X- a + e)- Vxx{X- a, 6)| = 0{e^'-^) . 



Due the above estimates it is easily seen that 7^ G r(A), Ve^. Due to 
Lemma [T7] and Corollary [H to prove the thesis it is sufficient to show that 
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the functionals Ki{\] a + e^) and K2{\] a + e^) converge in norm to i^i(A; a) 
and K2{\] a). Here Ki{\; a), i = 1, 2 are defined as in Corollary[T] We notice 
that the norm of the functionals are just the L^{'^i) norm of their integral 
kernels. 

We first consider the functionals A'2(A;a + e^). Due to the above esti- 
mates 

7 7 

A2Q;(/i, — )• A2a(/x), uniformly on 7i([0, 1]) as n — )• oo . 

Hence the sequence a(//, e^) converges in norm L^(7i) to and the 

sequence K2{\]e\) converges in operator norm to i^2(A;a). 
We consider now the sequence Ki{\\ a + e^). 

To prove the convergence of the above sequence of operators, it is suffi- 
cient to prove that 

gS,(Aia+£^)-Sfc(/.;a+4) ^ gSfc(A;0)-Sfe(M;0) uniformly on 7i([0, 1]) as ?i ^ oo . 
We first note that 



gSfc(A;0)-5fc(/i,;0) _gSfc(A;a+ej,)-Sfe(/.;a+ej,) ^ g5fc(A;0)-5fe(/.;0) ^ _ ^g{li;e)^ 

Using estimate (|42p . it is easy to show that g(fi;e) = f{e)0{fi^), where 
/(e) — 7- as e —7- and < 5 << 1. Therefore the difference of the expo- 
nential functions converges uniformly to 0. 

□ 

We can prove Theorem |2] (ii) and (iii). 

Indeed from (|29p . it is easily seen that (choosing one of the two branches 
of the gauge transform) 



M/fc(A;z) = G(A,z)-i$fe(A;z) 

if |argA- ^1 < ^ and |^ » 0. 
Moreover from (1291). it follows that 



lmi(z - a)$f ^(A, z) = iV2^[^\X; a) . 
Hence Theorem [2] (ii) and (iii) follow from Theorem [TUJ 
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